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Solid State Physics 


(Ch. 3) 


Last weeks: 

• Diffraction from crystals 

• Scattering factors and selection rules for diffraction 


Lecture 4 



Today: 

• Lattice vibrations: Thermal, acoustic, and optical 
properties 

This Week: 

• Start with crystal lattice vibrations. 

• Elastic constants. Elastic waves. 

• Simple model of lattice vibrations - linear atomic chain 

• HW1 and HW2 discussion 
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Material to be included in the 1 st QZ 

• Crystalline structures. Diamond structure. Pa 
7 crystal systems and 14 Bravais lattices 




d 


Crystallographic directions 

and Miller indices d m = 


Definition of reciprocal lattice vectors: 


'h 2 


7 2 A 


K a 2+ b 2+ c 2 ; 


1/2 


_ 2 71 III .2 71 

b, = — -a, xa., b. =-—a.xa, * b, = —a, xa, 

I y 1 “ Y - ' V 


• What is Brillouin zone 

• Bragg formula: 2d-sin 6 = mk \ Ak = G 
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•Factors affecting the diffraction amplitude: 

Atomic scattering factor (form factor): r _ f n ( r )e iAk ' r 'd 3 r 

reflects distribution of electronic cloud. a J 

In case of spherical distribution / = {4 nr 2 n(r) ' r ) 


A k-r 


dr 


•Structure factor 


-w —r ^ /* 2 tci (hu j kv j Iw j ) 

F = LLj e 

j 


where the summation is over all atoms in unit cell 



Na® •Cl 


• Be able to obtain scattering wave vector or frequency from 
geometry and data for incident beam (x-rays, neutrons or light) 
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Material to be included in the 2 nd QZ 
TBD 


Elastic stiffness and compliance. Strain and stress: definitions and 
relation between them in a linear regime (Hooke's law): 


/,= Z c » 


• Elastic wave equation: 




■ d U ^ eff 5 u 


“I®, 


ijkl ^kl 


kl 


dt 2 


p dx 2 


sound velocity v = 


C 


eff 


P 
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• Lattice vibrations: acoustic and optical branches 
In three-dimensional lattice with s atoms per unit cell there are 
3s phonon branches: 3 acoustic, 3s - 3 optical 


• Phonon - the quantum of lattice vibration. 

Energy hco\ momentum f )q 

• Concept of the phonon density of states 

• Einstein and Debye models for lattice heat capacity. 


Debye temperature 


ft 

e °'T. 


6ttN 

V 


, 1/3 


Low and high temperatures limits of Debye and Einstein models 


Formula for thermal conductivity 


K = -Cvl 
3 


• Be able to obtain scattering wave vector or frequency from 

geometry and data for incident beam (x-rays, neutrons or light) 
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Crystal Vibrations 


• The study of lattice vibrations 
is important for 
understanding how energy is 
absorbed in solids. 

• At low temperatures, most of 
the energy that is absorbed 
is of two types: lattice 
energy, and electronic 
energy - the heat can go to 
the electrons, or to vibrations 
of the positive charges in the 
lattice. 

• We will eventually show 

that 

this can be measured with 
the heat capacity: C = dU/dT 

• For a metal, C = yT + 6T 3 

Electronic term Lattice term 



C/T = y + PT 2 


U 0 



Elastic properties 

Elastic properties are determined by forces 
acting on atoms when they are displaced from 
the equilibrium positions 

R 


Taylor series expansion of the energy near 
the minimum (equilibrium position): 


U(R) = U 0 + 


dU 


dR 


(R-R 0 ) + 


1 d 2 U 


2 dR 2 


(R-R 0 ) + ... 


«0 


dU 


For small displacements, neglect higher terms. At equilibrium, 

OK 

S °’ U(R) = U 0 +^ ' 8U 


= 0 


where k = ■ 


dR 2 


u = R - R 0 - displacement of an atom from equilibrium position 
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dU 

force F acting on an atom: F =-= -ku 

SR 

k - interatomic force constant. This is Hooke's law in simplest form. 

Valid only for small displacements. Characterizes a linear region in which the 
restoring force is linear with respect to the displacement of atoms. 

Elastic properties are described by considering a crystal as a homogeneous 
continuum medium rather than a periodic array of atoms 

In a general case the problem is formulated as follows: 

• Applied forces are described in terms of stress cr, 

• Displacements of atoms are described in terms of strain s. 

• Elastic constants C relate stress crand strain e, so that cr = Cs. 

In a general case of a 3D crystal the stress and the strain are tensors 
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Stress has the meaning of local applied “pressure”. 

Applied force F (F x , F y FJ —> Stress components c^ ( i,j = 1, 2, 3) 

x = l,y = 2, z = 3 

General case for stress: i.e cr , 7 


c<j. 
F. =— B- 

1 °' x i 


■ dcr, 


\kdV = \-fdV = \da l dS J 


v dx j 



Shear forces must come in pairs: <jy = cr /( (no angular acceleration) 

=> stress tensor is diagonal, generally has 6 components 

10 


5 













Stress has the meaning of local applied “pressure”. 


Applied force F (F x , F y FJ —> 
General case for stress: i.e 


Stress components (ij = 1, 2, 3) 

a x = l,y = 2, z = 3 

Oj J-^L 

J dx : 


Hydrostatic pressure - stress tensor is equivalent to a 
scalar: i.e cr xx =c7 yy =cr zz 




-p 0 0 

0 -p 0 
0 0 — p 


Stress tensor is a “field tensor” that can have any symmetry not 
related to the crystal symmetry. Stress tensor can change the 
crystal symmetry 
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Stress has the meaning of local applied “pressure”. 

Applied force F (F x , F y FJ —» Stress components (/',/' = 1, 2, 3) 

x = 1, y = 2, z = 3 

Compression stress : i = /, i.e cr xx , cr yy , cr zz 


\| 

N 
\ 

\ 

V 

\ area A x 


Fx 


u 


<J rr — ■ 


A.. 


Shear stress : i +j, i.e <r xy , o yx , cr xz o zx , cr yz , o zy 


N 


O xx 


\ 

\ 


t F- 


CF — —- 

** A. 


I 




area A x +— 

Shear forces must come in pairs: cr :j = o ; 7 (no angular acceleration) 
=> stress tensor is diagonal, generally has 6 components 12 
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Strain tensor components are defined as 


£ ij 


Su 

dxj 


Compression strain e 


du x 

dx 


Shear strain ( e yx e xz e zx e yz ): 




<K 

dy 



Since o■■ and cr- always applied together, we can 
define shear strains symmetrically: 


£..=£..= — 

V J> 2 


du du- 

' + 1 


dxj fix 


So, the strain tensor is also diagonal and has 6 components 
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Elastic stiffness (C) and compliance (S) constants 

relate the strain and the stress in a linear fashion: 


This is a general form of the Hooke’s law. 




tjki £ ki 


kl 


£ U 


=2X 


kl 


6 components 6 £ tj —> 36 elastic constants 
Notations: C mn where 1 = xx, 2 = yy, 3 = zz, 4 = yz, 5 = zx, 6 = xy 
For example, C n = C xxxx , C n = C xxyy , C 44 = C yzyz 

Therefore, the general form of the Hooke’s law is given by 


ijkl&kl 


compression 


mixed 




fc,, 

Cn 

Cn 

Cn 

Cn 

c i« > 



°yy 


C 21 

C 22 

C 2 3 

C 2 4 

Q 5 

c 6 

£ yy 

compression 



Q 

C 32 

C33 

C34 

c„ 

Q, 

£ zz 




C n 

C 42 

C 43 

^44 

C 4 , 


£ r- 


<Tx 


C 51 

C 5 2 

C53 

C 5 4 

c 55 



shear 



Iq. 

c 62 


C 64 

C a 

^66 J 

K. € *y > 


stress 


mixed 



shear 


strain 
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Elastic constants in cubic crystals 

Due to the symmetry (x, y, and z axes are equivalent) c„ C 22 C 33 ; 

Q 2 = C 2J = C 13 = C 31 = C 23 = C 32 ; C 44 = C 55 = C 66 

Also, the off diagonal shear components are zero: 

C45 = C54 = c 46 = C 64 = C 56 = C 65 

and mixed compression/shear coupling does not occur: 

C45 = C54 = C 46 = C 64 = C56 = C 65 


=> the cubic elastic stiffness 

fc.. 

Ca 

C n 

C u 

Ca 

C n 


0 

\ 

tensor has the form: 




C 12 

C n 

C n 

C 44 

0 

0 

only 3 independent constants 


0 


0 

C„ 

0 





0 

0 

c 

'~44/ 























Elastic constants in cubic crystals 


Longitudinal compression „ _ cr xx _ F/A 

^11 — — 7 

(Young’s modulus): s xx u/L 



Transverse expansion: 


C, 


12 




F 


Shear modulus: 


C 


44 



F/A 

~T 
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Uniaxial pressure setup 

for optical characterization of correlated oxides 


Pressure control 


Optical 

cryostat 
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Elastic waves 

Considering lattice vibrations three major approximations are made: 

• atomic displacements are small: u« a , where a is a lattice 
parameter 

• forces acting on atoms are assumed to be harmonic, i.e. 
proportional to the displacements: F=- Cu 

(same approximation used to describe a harmonic oscillator) 

• adiabatic approximation is valid - electrons follow atoms, so that 
the nature of bond is not affected by vibrations 

The discreteness of the lattice must be taken into account 
For long waves A » a, one may disregard the atomic nature - 
solid is treated as a continuous medium. 

Such vibrations are referred to as elastic (or acoustic) waves. 

19 


Elastic waves 

First, consider a longitudinal wave of compression/expansion 

a 

mass density p 
segment of width dx at 
the point x; 

elastic displacement u 


ZlZ! 



1 

— » 

1 

F(.y) 

1 

“ 

* 

F(.v+d\ ) 


x .v+dv 


777 ^ - =Yj F -> (7 0Adx) ~~~t = F(X + dx) - F(x) 

dt dr 


d 2 u 1 dF Sex 

=> p —- =-= —“ 

dt 2 A dx dx 


where F/A = o YY c -£=.=££i 

XX J L — - /r 



Assuming that the wave propagates along the G 

[100] direction, can write the Flooke’s law in the form xx 

du 

Since £ xx = get wave equation: 

dx 


d 2 u _ C n d 2 u x 
dt 2 p dx 2 


F 
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A solution of the wave equation - longitudinal plane wave 

u{ X, t ) = where q - wave vector; frequency co = v L q 


v 


L 



- longitudinal sound velocity 


Now consider a transverse wave which is controlled by shear stress 
and strain: « 


In this case 




cfu _ da v 
^ dt 2 dx 

=> wave equation is 


x -v+dx 

where cr AV = C 44 s xy and e xy 


d 2 u _ C 44 d 2 u x 
dt 2 p dx 2 



_ du 
dx 

- transverse 
sound velocity 
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Two independent transverse modes: displacements along y and z 
For q in the [100] direction in cubic crystals, by symmetry the 
velocities of these modes are the same - modes are degenerate 

Normally Q/ > C 44 => 'T > v T 

We considered wave along [100]. 

In other directions, the sound velocity depends v = 
on combinations of elastic constants: 

C eff - an effective elastic constant. For cubic crystals: 


Mode 

q || [100] 

q II [lio]. 

q II [HI] 

L 

Cn 

1 /2(Cn+(?i2+2C44) 

1 /3(Cn+2 Ci 2+4C 44 ) 

Ti 

C 44 

C44 

’/3 (Cn—C12+C44) 

t 2 

C 44 

V2(C n -Cn) 

V3 (Cn—C12+C44) 


Relation between co and q - dispersion relation. For sound co =vq 

22 
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Model of lattice vibrations 

one-dimensional lattice: linear chain of atoms 

U„-\ "n «»+1 

harmonic approximation: 

force acting on the n th atom is , , ◄-► 

;7-l 77 77+1 a 

F n =C(M„ + i-W„)+C(W„_i-«„) 

equation of motion (nearest neighbors interaction only): 

d 2 u 

M ~3t 2=Fn=C ~ Un) + C ~~ U n ) = ~ C ( 2U n ~ U n + 1 “ U n -1) 

M is the atomic mass, C - force constant 
Now look for a solution of the form u(x, t ) = Ae i( ' qXn ~ a *'* 

where x n is the equilibrium position of the n-th atom ~^x n = na 

obtain M{—Q}~')e iqtta = — — g'ft"- 1 )** J 23 


M{-co 2 )e iqna = -C[2^" < ’ -e ,4( ” +1)< ’ -e , ‘ 7(B ’ 1)o ] 


-> 


A/<y 2 =r(2-«^ -«*") = 2C(l-cos?«) = 4Csin 2 ^- 


_ 4C 

the dispersion relation is 0} ~x~ 


sin 


qa 


Note: we change q —>q + 2n/a the atomic displacements and 
frequency co do not change => these solutions are physically identical 


can consider only 
K 71 

- <q < — 

a a 


i.e. q within the first Brillouin zone 3 
The maximum frequency is 2^/^ 0 0 
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Dispersion relation 

• Another way of showing this, physically, is by looking at the 
wavelengths of vibration: 



• The wave represented by the solid curve conveys no new 

information not given by the dashed curve. Only wavelengths longer 
than 2a are needed to represent the motion. This means that K is 
limited by -TT/a < K < n/a 
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Phase and group velocity 


co 

phase velocity is defined as v p = — 


dco 

group velocity v = — 

dq 



v g = 0 at the boundaries of the Brillouin zone (q = ±nId) => 
no energy transfer - standing wave 
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Features of the Dispersion Relation 


Most of the waves 
described by the 
wavevedor K are 
travelling waves (meaning 
that the waves propagate 
through the lattice). 
However, at the zone 
boundary, which is 
K = +/- n/a the wave 
becomes a standing 
wave . This means that 
the wave itself moves 
neither to the left or two 
the right. 

If is also a standing wave 
at K = 0 (the wavelength 
is infinite ) 





(at Ka = +/-tt these are vibrating out of phase, and so the wave 
does not propagate) 


27 
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Diatomic lattice 


one-dimensional linear chain, atoms of two types: M ; and M 2 

Mi Mi 

■4 -► 

n —1 n n +1 a 


OJ 



Optical Phonons 
can interact with light 

For diamond 

Optical phonon 
frequency is * 1300 cnv 1 

X» 7700 nm 
(far-IR) 
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Model of diatomic lattice 

one-dimensional linear chain, atoms of two types: M } and M 2 

Mi M 2 

Treat in similar way, but 

need two equations of ^-► 

motion: n W+1 a 

M \ = ~ c ( 2w » “ w » + i -*Vi) 

a/ 

w 2 i 3 lL= - C(2u ”*'-"^-“.) 

at 


Again, look for a solution of the form U n — A x e 


— A J(qna-cot) 


_ a i(q(n+\)a-cot ) 


U n + 1 = A 2 e 


30 


Substitute this solution into equations of motion 
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get system of two linear equations for the unknowns A 1 and A 2 
In matrix form: 

/i — 11/1 tfI —/a i’ll 's.ii/i n. 


"2 C-M,co 2 

-2 C cos qa 

~A 

-2 C cos qa 

2C -M^a) 2 

A. 


determinant of the matrix must be zero 


(2C - M l or )(2C- M 2 co 2 ) - 4C 2 cos 2 qa = 0 
Solve this quadratic equation, get dispersion relation: 


' i 

h l ] 

±c. 

V 

l 

1 1 
l__ 

2 

4sin 2 qa 

A,' 


a' 

m 2 J 



Depending on sign in this formula there are two different 
solutions corresponding to two different dispersion curves 



q 

The lower curve - acoustic branch, the upper curve - optical branch. 
at q = 0 for acoustic branch co 0 = 0; A 1 = A 2 
=> the two atoms in the cell have the same amplitude and the phase 
dispersion is linear for small q 


for optical branch 

at q = 0 




1 J_ 

M, M 


2 J 


MjA i +M 2 A 2 — 0 


=> the center of mass of the atoms remains fixed. The two atoms move out of 
phase. Frequency is in infrared - that's why called optical 32 
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Transverse optical and 
acoustical modes 



acoustical mode 


• The acoustic branch has this name because it gives rise to long wavelength 
vibrations - speed of sound 

• The optical branch is a higher energy vibration (the frequency is higher, and 
you need a certain amount of energy to excite this mode). The term “optical” 
comes from how these were discovered - notice that if atom 1 is +ve and 
atom 2 is -ve, that the charges are moving in opposite directions. You can 
excite these modes with electromagnetic radiation (ie. The oscillating electric 
fields generated by EM radiation) 


Acoustic/optical branch 



34 
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Zone boundary 

• At the zone boundary, Ka = +/- tt, so cos(Ka) = -1 

• Substitute this back into our original equation: 

„ 2 2 C(M l+ M 2 ) ^ V( 2C Wi + M 2» 2 -4M 1 M 2 (-4C 2 ) 

UJ — 

2 M,Mj 

1 Miir *£)]“ 

Roots: 2 2C 

CO =- 

M 1 

2^ 2C 

CO =- 

m 2 

2 M x M 2 

Optical phonon branch 

n 

^ (2 C/M 2 ) W2 

Mi > M 2 

_ _ C/M,) 1 * 

X 

Acoustical 
yX phonon branch 

x 

X i 


If - K 

a 

35 


Two Atoms/Primitive Basis 


* If there are p atoms in the unit 
cell, there has to be 3p branches 
in the dispersion relation: 3 
acoustical branches, and 3p-3 
optical branches (there are more 
ways for the atoms to move for 
the optical modes) 

* So T for Germanium, which has 2 
atoms/unit cell (diamond 
structure: 000 and !4 % % ), 
there has to be 6 branches: one 
longitudinal acoustic (LA), 2 
transverse acoustic (TA) h one 
longitudinal optical (LO), 2 
transverse optical (TO) 



(phonon dispersion curve in the [111] 
direction for Germanium at 80 K, as 
determined by inelastic neutron scattering) 

36 
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• How do we measure dispersion curves? Inelastic neutron scattering 

• If you can measure the energy lost by the neutron (by causing a vibration in the 
solid), and you can measure which direction you created the wave (the 
wavevector), then you can construct a dispersion curve 
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Summary 

❖ Elastic properties - crystal is considered as continuous anisotropic 
medium 

❖ Elastic stiffness and compliance tensors relate the strain and the 
stress in a linear region (small displacements, harmonic potential) 


Hooke's law: 


°ij 


C ijkl £ kl 


£ ij 


2X 


iju^u 


kl 


kl 


♦> Elastic waves q 2 u C eff d 2 u sound velocity 


v - 


if 


c. 


eff 


P 


dt p dx 

❖ Model of one-dimensional lattice: linear chain of atoms 

❖ More than one atom in a unit cell - acoustic and optical branches 

❖ All crystal vibrational waves can be described by wave vectors 
within the first Brillouin zone in reciprocal space 


What do we need? 3D case consideration 

Phonons. Density of states 


38 
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HW1 


Compute the theoretical density of NaCl based on its crystal structure. 



For NaCl structure, the 
crystal lattice parameter is 
a= 2 ( r Na + + rCI -), 
where r is ionic radius. 


_ jM _ 4(A Na + A C1 ) 
P "V a‘N A 
4ions (22.99 + 35.45) g/mol 


2 [(0.102x 10‘ 7 + 0.181 x 1 O' 7 )] cm 3 (6.023x 10 23 )ions/mol 


= 2.14g/cm 3 ■* - 

(actual = 2.16g/cm 3 ) 


39 
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Vibrations in three-dimensional lattice. 


Phonons 

Phonon Density of states 
Specific heat 



(Ch. 3.3-3.9) 
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Three-dimensional lattice 

In simplest ID case with only nearest-neighbor interactions we had 

equation of motion 

solution 

d 2 u 

M-^ = F n =C(u n+l -u n ) + C(u n _ 

u(x,t) = Ae' iv --“ ) 

In general 3D case the equations of motion are: 

/ / Atom a 

m 5 Un “ = y F m/J 

a o,2 / , na 

Ot in,/? 

j/- \ 

A/ unit cells, s atoms in each > 

f nth elementary cell 

3/V’s equations 

\ i n = (n 1 ,n 2 ,n 3 ) 

Fortunately, have 3D periodicity => 
Forces depend only on difference m-n 

\ I fna = + r„ 

Write displacements as 

r n = n,a, + + n 3 a 3 \ I 

v M « 

* (0,0,0) 


substitute into equation of motion, get 


-ryV 


od-IZ 






=F m ^ J e iq '^ m 


-**„) 


u, 


Mr 


,(q) = o 


Dj' (q) -dynamicalmatrix 


- ® 2 w«, (q) + X Dg (q )u g (q) = 0 

P,j 

=> Det{z)f(q)-« 2 l}=0 

- dispersion relation 

3s solutions - dispersion branches 

O 

3 acoustic, 3s - 3 optical 

direction of u determines polarization 
(longitudinal, transverse or mixed) 

Can be degenerate because of symmetry 


phonon dispersion 
curves in Ge 
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Phonons 

• Quantum mechanics: energy levels of the harmonic oscillator are quantized 

• Similarly the energy levels of lattice vibrations are quantized. 

• The quantum of vibration is called a phonon 

(in analogy with the photon - the quantum of the electromagnetic wave) 

Allowed energy levels of the harmonic oscillator: E = (ll + 'A)hO) 

where n is the quantum number 

A normal vibration mode of frequency co is given by u = 

mode is occupied by n phonons of energy ha>; momentum p = hq 

Number of phonons is given by Planck function: ^ 

(T- temperature) n = f - 

ghw/kT _ 1 

The total vibrational energy of the crystal is the sum of the energies of the 
individual phonons: 

(p denotes particular E = ^E V , = 

phonon branch) w w 43 


Phonon Momentum 


• You might think that these phonons have momentum which 
is p = hK, just like photons have momentum p = E/c. 

• Be careful! Phonons don’t carry momentum like photons 
due. They can interact with particles like they have a 
momentum (for example, a neutron can hit a crystal, and 
start a wave by transferring momentum to the lattice). 

• However, you have to think of this momentum as being 
transferred to the whole lattice . The atoms themselves are 
not being translated permanently from their equilibrium 
positions. 

• The only exception to this rule is the K = 0 mode, where the 
whole lattice translates. This, of course, carries momentum. 
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Phonon Momentum 


• For all practical purposes, a phonon acts as if it carries a 
momentum p = hK, which is sometimes called the crystal 
momentum. 


We have already seen examples of where this is used. 
For example, if an x-ray interacts with a lattice, then we 
know that the scattered ray(k’) and incident ray (k) have 
to be related by: 


k' = k + G 


Where G is a reciprocal lattice vector. In this process, 
the whole crystal recoils with momentum -hG 



Crystal 


Recoil -G 
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Phonon Momentum 

If this photon interacts inelastically with the lattice, then 
we have: 

(phonon is created) 


k’ + K = k + G 


Or: 


k’ = k + K + G 


(phonon is absorbed) 



K phonon absorbed 
"k 



Crystal 
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Density of states 

Consider ID longitudinal waves. Atomic displacements are given by: 

u = Ae iqx 

Boundary conditions: external constraints applied to the ends 
Periodic boundary condition : u(x = 0) = u(x = L) 


Then e iqL = 1 => condition on the admissible values of q\ 

2n 

q = —n where n - 0, ±1, ±2,... 

L 


regularly spaced points, spacing 2n/L 

Number of modes in the —dq 
interval dq in g-space : 2 n 

Number of modes in the frequency 
range (Q, o + rfco): D( )d(0 = _L d 

In 

D( co) - density of states 
determined by dispersion co = co(q) 



D(co) - ——— 
n dco /dq 


Density of states in 3D case 

Now have u = \ e ^ qxX+q > y+q - ^ 

Periodic boundary condition : e lc,xL = e qyL = e' q:L = 1 


, \ f,2n In 2n^ 

3 T' m T’"T 

Plot these values in a g-space, 
obtain a 3D cubic mesh 

number of modes in the spherical 
shell between the radii q and q + dq\ 


\2 n j 


4nq 2 dq = 


V 


( 2 n) 


■4nq 2 dq 


V = L 3 - volume of the sample 

Vq 2 

Density of states D(co) = 


l,m,n- integers 


02 contour 


2 n 2 dco/dq 



dq 


24 




















Few notes: 

• Equation we obtained is valid only for an isotropic solid, 

(vibrational frequency does not depend on the direction of q) 

•We have associated a single mode with each value of q. 

This is not quite true for the 3D case: for each q there are 3 different 
modes, one longitudinal and two transverse. 

• In the case of lattice with basis the number of modes is 3s, 
where s is the number of non-equivalent atoms. 

They have different dispersion relations. This should be taken into 
account by index p =1 ...3s in the density of states. 


Lattice specific heat (heat capacity) 

dO (dF ^ 

Defined as (per mole) C = — If constant volume V C v = 

dT v 


\OT Jy 


The total energy of the phonons at temperature Tina crystal: 


£ _ y l n /g'j _ o (the zero-point energy is chosen as the 

“ ' qp / p origin of the energy). 


1 


n) = 


fico/kT 


-1 


Z ha>J q) 

—;- - 

ha> (q)/k B T 


qp 


-1 


replace the summation over q by an integral over frequency: 

E^dcoD^a,)- hC ° 


tico/koT i 

e 8 -1 


Then the lattice heat capacity is: 


Cr- 


tiro 

k&r j 


tia/koT 


8T 

Central problem is to find the density of states 


sftcoik B T _j^2 
50 


25 













Debye model 


•assumes that the acoustic modes give the dominant contribution 
to the heat capacity 

•Within the Debye approximation the velocity of sound is taken a 
constant independent of polarization (as in a classical elastic 
continuum) 


The dispersion relation: co =vq, v is the velocity of sound. 
In this approximation the density of states is given by: 


D(g>) = 


Vq 2 


1 


Vq 2 1 Vco 2 


In 2 dcoldq 2 n 2 v 


2k 2 v 3 


Need to know the limits of integration over co. The lower limit is 0. 

How about the upper limit ? Assume N unit cells is the crystal, only 

one atom in per cell => the total number of phonon modes is 3N => 

®d f r 2 3 Ar y / 3 Debye 

6 7T V N /, 2 V/3 

co n = - = v(6 n n J 


IJ D{co)d(o = 3 N 


p o 


V 


frequency 


The cutoff wave vector which corresponds to this frequency is 


q D = — = 

v 


( 6tt 2 N Y /3 m °des of wave vector larger than q D are not 

allowed - number of modes with q <q D 
exhausts the number of degrees of freedom 


V 


Then the thermal energy is 
where is "3" from ? 

3 Vh C0D 


C0 D T/ .,2 

E = 3 f dco-^—r- 

{ 2tt 2 \ 


tuo 


3 gfio/kgT _ | 


„ 3 Vh 7 , ax’ 3 VkiT 4 x ? , x 3 

E = —aco f _ -= — , , , ax - 

2^ 2 v 3 l e t,a/k » T -1 27T 2 v 3 h 3 l e x -1 


where x = hco/k B T and x D = hco^kgT = Oj/T 


Debye temperature: 


a -!H 

D k 
k b v 




67tN 
V 


1/3 
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The total phonon energy is then E = 9 Nk B T 


7 - \ X D 3 

- 

0 n i e x -l 


T_ 

y° D j 0 

where N is the number of atoms in the crystal and x D = 0 L /T 
To find heat capacity, differentiate 


J .r , CO D 3 

r f j C0 

E =-- —- a CO —;— rr-= — 

^ 3 I no* k n T i 

2 n v i e 5 -1 


So, 


cv = 


3 Vh f co e * 

-—-r dco -- = 9Nk R 

InWkJ 1 { (e halk ' T -\) 


2 C0 D 


4 Tico/knT 

co e 


’~r \ D 

- f* 

\ e o) l 


4 x 

x e 


(*'->) 


In the limit T»6 D , x « i,=> C v = 3M 5 - Dulong-Petit law 
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Einstein model 

The density of states is approximated by a ^-function at some co E : 
D(E) = NS(co - co E ) where N is the total number of atoms - 

simple model for optical phonons 

„ 3 Nti<o F 

Then the thermal energy is E = — — -r ^— 

he> E lk„T _ i 


The heat capacity is then 


C y =—=3M, 

cT 


V J 


/ k a T 


^ e hco E /k B r 


The high temperature limit is the same as that for the Debye model: 
C v = 3 Nk B - the Dulong-Petit law 


At low temperatures C v ~ e~ hm/k “ T - different from Debye V law 

Reason: at low T acoustic phonons are much more populated => the 
Debye model is much better approximation that the Einstein mocfel 


Real density of vibrational states is much more complicated than those 
described by the Debye and Einstein models. 

This density of states must be taken into account in order to obtain 
quantitative description of experimental data. 


The density of states for Cu. 

The dashed line is the Debye 
approximation. 

The Einstein approximation 
would give a delta peak at 
some frequency. 





oi, 10 13 rad/s 


28 














Summary 

❖ In three-dimensional lattice with s atoms per unit cell there are 
3s phonon branches: 3 acoustic, 3s - 3 optical 

❖ Phonon - the quantum of lattice vibration. 

Energy hco\ momentum Fiq 

❖ Density of states is important characteristic of lattice vibrations; 

It is related to the dispersion co = cofq). 

Vq 1 

Simplest case of isotropic solid, for one branch: D(co) = --— 

2n dcojdq 

❖ Heat capacity is related to the density of states. 

❖ Debye model - good when acoustic phonon contribution dominates. 
At low temperatures gives C v oc T 3 

❖ Einstein model - simple model for optical phonons (co(q) is constant) 

❖ At high T both models lead to the Dulong-Petit law: C v = 3Nk B 

❖ Real density of vibrational states is more complicated 57 




